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Recent experiments showing scaling of the intrachromosomal contact probability, P{s) ^ s~^ with 
the genomic distance s, are interpreted to mean a self-similar fractal-like chromosome organization. 
However, scaling of P{s) varies across organisms, requiring an explanation. We illustrate dynamical 
arrest in a highly confined space as a discriminating marker for genome organization, by modeling 
chromosome inside a nucleus as a homopolymer confined to a sphere of varying sizes. Brownian 
dynamics simulations show that the chain dynamics slows down as the polymer volume fraction 
(0) inside the confinement approaches a critical value <fc- The universal value of ~ 0.44 for 
a sufficiently long polymer {N ^ 1) allows us to discuss genome dynamics using 0 as a single 
parameter. Our study shows that the onset of glassy dynamics is the reason for the segregated 
chromosome organization in human [N 3 x 10^, (f> > 4>T)^ whereas chromosomes of budding yeast 
{N ^ 10®, (f < 0^) are equilibrated with no clear signature of such organization. 


Chromosomes exhibit dramatic changes in their 
spatial organization along the cell cycle. In the 
metaphase, they are condensed into compact blob-like 
structures [T], whereas in the interphase they decon- 
dense to a less compact coil-like structures. Interphase 
chromosomes are not random but form territories [2], 
and their organization may be fractal-like [3]. Ad¬ 
vances in experimental techniques Hd have provided 
quantitative details of chromosome organization in the 
form of chromosomal contact maps describing how dis¬ 
tant loci are structurally organized. The contact prob¬ 
ability of two loci separated by a genomic distance s 
scales as P{s) ^ differing from P{s) ^ in 

equilibrated polymer melts. The deviation of the ex¬ 
ponent from —1.5 is taken as an evidence that chro¬ 
mosomes form a non-equilibrium globule with segre¬ 
gated domains rather than a fully equilibrated globule 
with entanglements EE]. Such an interpretation of 
the structural organization based solely on P{s) is not 
universally accepted [9]. In addition, genome struc¬ 
ture could vary depending on the extent of maturity 
of human cells HO]. Still, the scaling of P{s) varies 
depending on organisms. It is therefore important to 
develop a theoretical framework for distinguishing be¬ 
tween genome structure in different organisms. 

From a biological perspective, it could be argued 
that the hierarchical and scale-free organization of 
chromosome, without knots, is beneficial for access 
to a target locus m or for the faster response to 
an environmental change by easing the condensation- 
decondensation process [T2j[T3]. Although the origin of 
chromosomal territories is controversial because equi¬ 
librium polymer configurations with many loops nat¬ 
urally produce segregated domains as well [HI [15] , a 


major non-equilibrium effect, glassy dynamics of the 
genome under strong confinement, should not be over¬ 
looked from a contributing factor in chromosome fold¬ 
ing. The relaxation time of a polymer via disentan¬ 
glement [16] {Trep ^ [121 [ini Ull HE]) could be 

far longer, effectively permanent for higher organisms, 
than the cell cycle time {tccIi) [HI US] for a large N. 
Furthermore, a substantial increase of polymer relax¬ 
ation time is also expected in a strong confinement as 
is the case for DNA inside viral capsid m even when 
N is not too large. Thus, to fully describe the genome 
structuring, it is imperative to understand the poly¬ 
mer dynamics under confinement and how it might 
vary across various species. The major goal of this 
work is to develop a physical basis, using relaxation 
dynamics as a quantitative measure, to discriminate 
between genome organization in different organisms. 

Although explicit models that consider circular 
DNA or multi-chains and specific contacts based on 
Hi-C contact maps [2QH22] are possible, here we study 
the dynamics of homopolymers confined to a sphere 
of varying sizes as a first step towards understand¬ 
ing the dynamical features of interphase chromosomes. 
We consider a single self-avoiding polymer chain, rep¬ 
resenting chromatin fiber, confined to a sphere and 
employ dynamical measures previously used to study 
supercooled liquids [23ll25] as a vehicle to investigate 
non-equilibrium effects. Our major finding is that the 
dynamics and organization of homopolymer vary dra¬ 
matically as the extent of confinement is increased. 
When this result is translated into genome organi¬ 
zation, we find that bacteria and yeast chromosome 
folding can be thought of as an equilibrium process 
whereas glassy behavior governs the territorial organi- 
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FIG. 1: Polymer {N = 300) confined to spheres, (a) Snapshots from simulations. The value of potential energy (scale on 
the left) for each monomer shows that the spatial heterogeneity increases as (j) approaches 0c(300) ~ 0.404 (see Fig.S5 
for other snapshots). Contact maps from three distinct polymer conhgurations at 0 = 0.402 near 0c(300) are shown on 
the right, (b) Flory laws Rg ~ and ~ are satisfied for unconfined and strongly confined (0 = 0.402) chains 

of varying N, respectively, (c) RDFs (see SM) at three 0 values, (d) Inter-segmental contact probabilities (see SM for 
dehnition) with increasing 0 from 0 to 0.404 (N = 300). (e) Distributions of monomer energy (e^ = 

where the surface interaction term is excluded from the calculation. See SM), P(e) = N~'^ ~ increasing 

0. Divergence of the standard deviation (cr^ = (e^) — (e)^ where {e^) — e^P{e)de) near 0c is shown in the inset. 


zation in humans. These inferences cannot be drawn 
from genome contact maps alone, which has been the 
sole focus on chromosome folding. 

The equilibrium aspects of confined polymers are 
well understood m [26]. The equilibrium free en¬ 
ergy of polymer confined to a sphere is not exten¬ 
sive [23 EH] in contrast to polymer localization in a 
slit or a cylinder. Furthermore, as the extent of con¬ 
finement increases, the volume fraction, defined by 
cf> = {Rl/Rsf (Fig[^, see Supplemental Material) in¬ 
creases, and more importantly, the equilibration time 
of the chain (reg) increases dramatically. If r^q for a 
genome is longer than finite cell doubling time {tccIi) 
then the decondensation-condensation cycle dynamics 
of the genome should be under kinetic control. We 
explore these aspects in the context of genome fold¬ 
ing using simulations of homopolymers confined in a 
sphere (see SM for details), highlighting the confine¬ 
ment effect on polymer leading to the ultraslow glassy 
dynamics, such that Teq ^ Tceih which we will show is 
the case in human chromosomes {N ^ 10^) and viral 
DNA {N ^ 10^). 

In general, it is difficult to distinguish between non¬ 
equilibrium conformation of a polymer from its equilib¬ 
rium counterpart because polymer configurations for 


both cases could be similar, just as is the case for liq¬ 
uids and glasses. Indeed, the polymer size with in¬ 
creasing N satisfies the Flory relationship, Rg ^ 
with ly ^ 3/5 and 1/3 for weak and strong confine¬ 
ment, respectively (Figj^), crossing over the regime 
Rg ^ at 0^^^ ~ 0.2 where repulsion due to ex¬ 

cluded volume is counter-balanced by the confinement 
pressure. Therefore, in strong confinement Rg scal¬ 
ing cannot distinguish between equilibrium and non¬ 
equilibrium globules. The radial distribution function 
(RDF) between monomers at high 0 is reminiscent of 
the closely packed structure (Fig{^), suggesting that 
extent of confinement controls the chain organization. 

The scaling exponent a of the contact probability 
between two sites separated by the chain contour s, 
P{s) ^ s is one way to assess the chain organiza¬ 
tion (Alternatively, the average distance between two 
loci separated by 5, R{s) ^ can be used um- 
See Fig.S3). P{s) s is expected for unconfined 
self-avoiding walk (SAW) (see SM) [29H3T]. For an 
equilibrium globule under strong confinement, poly¬ 
mer chains are in near ©-condition because of the ef¬ 
fective cancellation between attraction and repulsion. 
Hence, we expect that P{s) ^ s~^ with a = 1.5 
[I7l|32|. In the case of strong confinement, however, 
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P{s) ^ s~^ in the range of s/a ^ 0(10) for N = 300, 
similar to the scaling observed in the Hi-C analysis 
of the chromosome in interphase 13 ES]. The range 
of s“^-scaling increases as the extent of confinement 
increases (Fi gUt)- 

However, P{s) scaling is not an indicator of the un¬ 
derlying dynamics. Even for a SAW chain with no spe¬ 
cific attractive interaction, dynamics can be arrested 
in strong confinement, preventing a full equilibration 
of the chains on relevant time scales. We document 
the emergence of glassy behavior under strong confine¬ 
ment by first calculating the potential energy of each 
monomer (Figsj^, [^, S5). The spatial heterogeneity 
of monomer energies at 0 = 0.402 is striking (Figsj^, 
S5), which is also indicated by the abrupt changes in 
the monomer energy distribution P(e) (Figj^) and the 
standard deviation inset). 

In the absence of obvious symmetry breaking, it is 
useful to characterize the dynamics using van Hove 
correlation function to discern the onset of glass-like 
behavior [23]. The correlation function, 

1 ^ 

i=i 

provides dynamical information of how the system re¬ 
laxes from its initial configuration, where fj{t) is a 
position of monomer at time t. The ensemble- 
averaged isotropic self-intermediate scattering func¬ 
tion {Fq{t)) is estimated by integrating Fq{t) over 
space with q = \^ and ai q = g'max = 27r/rs, where Vg is 
the position of the first peak in the total pair distribu¬ 
tion function (see Figj^). The onset of the structural 
glass transition is described by the density-density cor¬ 
relation function {Fq{t)) as a natural order parameter, 
which decays to zero in the liquid phase, but saturates 
to a non-zero value in the glassy phase even at long 
times. Thus, {Fq^^^{t)) provides information of how 
rapidly the polymer confined to a sphere loses mem¬ 
ory of the initial configuration (FigJ^). From physical 
considerations, {Fq^^^{t)) should vanish at long times 
{t oc) for 0 < 0c; the decorrelation time of the 
polymer configuration increases sharply as the extent 
of confinement (or 0) approaches its dynamical arrest 
value. (Pgmax(t)) at various 0 is well fit by a stretched 
exponential function ^ ^ and the dependence 

of Tc^ on 0 for different N (Figj^ is analyzed using the 
relation, 

r,(0;7V)=ro(7V)(0c(7V)-0)-"^. (2) 

The relaxation time rQ,(0; A") increases with 0 and di¬ 
verges at 0c(A). The stretching exponent 0 decreases 
with 0 (Fig.S4), in consistent with our findings in 



FIG. 2: Polymer dynamics under confinement probed using 
Ta(0). (a) (Pgmax(^)) with Varying 0 for iV = 300. The time 
on the abscissa is scaled by r = P jD (see SM). (b) (top) 
For polymer with A, Ta(0; A) are fit to Eq[^ (bottom) To 
obtain the universal scaling exponent of Tq, near 0c(A), the 
fit was made using Ta/To(A) = (0c(A) — 0)~^^ for all A, 
which confers Vr — 0.65. (c) Finite size scaling to obtain 
0^ = 0c(A > 1). 0c(A)s fitted to 0c(A) = 0?° - aA-"' 
give 0^ = 0.449, a = 0.552, and A = 0.432. 

Fig{^ that the system becomes more glassy as 0 in¬ 
creases. The set of rc(0;A), for various A, are de¬ 
scribed by a universal curve, satisfying log(rcK/ro) = 
—Vr log (0c(A) — 0), and hence we obtain a universal 
scaling exponent Vr ~ 0.65 for the dynamical arrest. 
The critical volume fraction 0c (A) is A-dependent but 
saturates to a finite value 0^ in the limit A ^ oc. 
From finite size scaling (Figj^), we obtain 0^ « 0.45. 

As an alternative to rc(0), the fluctuations in 
Fgmax (^)5 namely the generalized susceptibility X 4 (t) 
corresponding to the variance in Fq^^^{t), can distin¬ 
guish between the states below and above 0c clearly. 
The fourth order dynamic susceptibility m, used to 
quantify dynamic heterogeneity in structural glasses, 
is given by 

X4W = N [(F,_(t)2) - (F,_(t))2] . (3) 

The amplitude of X 4 (t), x™^^? increases with 0 
(FigjSgi), and the divergence of near (j)c{N) 

can be described using x™^"^(0; A) = X4(^)(0c(A) — 
0)-^x. The scaling exponent for dynamical arrest 
transition is found to be v^ 0.37. In the A ^ oo 
limit, 0^ = 0.44 (FigjS^), which is consistent with 
the 0^ = 0.45 from the analysis based on Eq§ 

The significance of the key finding that 0^ ~ 0.44 
becomes transparent by predicting the consequences 
for chromosome dynamics in various organisms. With¬ 
out confinement or any special interactions mediated 
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FIG. 3: Polymer dynamics under confinement probed us¬ 
ing (a) with varying 0 for iV = 300. 

(b) (top) = xl{N){4>c{N) - 4>)~''^ diverges at 

iV-dependent critical volume fraction (l)c{N). (bottom) 
The data of are rescaled with tq{N) and the ht us¬ 
ing Xa^^/X%{^) — {<t>c{N) — for all N ^ which gives 

0.37. (c) Finite size scaling: (t)c{N) ht to (t)c{N) = 
—aN~'^ gives = 0.443, a = 0.571, and ly" = 0.453. 
Note that u ^ v" (Fig. 2c). 


by proteins, genome occupies a large volume V ^ 
^{Rg)^ with Rg ^ lp{N/g)^^^ {Ip ^ 50 nm = g x 0.34 
nm/bp, thus g ~ 147 bp). Given that the nuclear 
sizes are similar (^ 0{1) /im), there could be a large 
variation in the nuclear volume fraction for different 
organisms that have different genomic size, N. 

(i) For bacteria (A^ = 10^ bp), Rg ^ 10 gm is greater 
than the bacterial cell size ^ 1 gm. As 1 bp corre¬ 
sponds to 1 nm^ [34] , the volume fraction for bacterial 
genome is 0 bac = 1 nm^/bp x 10 ^ bp/l/im^ = 10 “^ ^ 
0 ^, implying that glassy effects are not relevant. 

(ii) In eukaryotes, DNA chains are organized in nu- 
cleosomes. Thus, it is more appropriate to estimate 
the volume of chromosomes in terms of the number of 
nucleosomes rather than volume of bare DNA. Since 
each nucleosome, whose volume is I4uc ~ tt x 10^ x 3 
nm^ (15-20 nm width, 2-3 nm height), is wrapped by 
^150 bp DNA with a 50 bp-spacer between the neigh¬ 
boring nucleosomes [35|, 200 bp-DNA is required to 
compose one nucleosome. For budding yeast N ^ 10^ 
[36] , the volume occupied by the entire nucleosomes is 

= 10^ bp/200 bp X Kuc - 0.47 gm^; and the 
yeast nucleus volume is ^ ^ gm^. Therefore, 

0yeast = ~ 0.12, which is Smaller than 

0^. This explains the intrachromosomal contact fre¬ 
quency P{s) for yeast genome, pointing to an 

equilibrium globule jS] [371439] . 

(iii) Human nucleus size varies depending on the cell 


type and the stage of development, which results in 
P{s) ^ for stem cell and P(s) ^ 5 for ma¬ 

ture cell [9]. For illustrative purposes, we adopt the 
nucleus volume ^60 — 110 gm^ from the av¬ 

erage size of mammalian cell nucleus 2 x R ^ b — 6 
gm laiiD]. Since the volume taken by the entire 46 
chromosomes, as a diploid with 2 x AT ^ 2 X 3 X 10^ bp, 
is = 6 X 10^ bp/200 bp x Kuc - 3 x 10^^ nm^ 

the volume fraction of human genome is (/)human = 

y^human/yhuman _ q 3 _ q 5 > 0 ^. Qf particular 

note is that Shuman ^ 0^- Thus, the lack of nuclear 
space in human cell makes the chromosome dynam¬ 
ics intrinsically glassy, indefinitely slowing down the 
relaxation of chromosome configuration. This crucial 
conclusion based on the simple estimate of Shuman sug¬ 
gests the decondensation-condensation process, driven 
by a panoply of partner enzymes, is likely to be under 
kinetic control. 

(iv) The volume fraction of DNA (L ^ 6 — 60 gm) 
inside a viral capsid {R ^ 25 — 50 nm) using L ^ 30 gm 
and P ^ 35 nm is ^vims ~ 0.5 > (l)c [T51[iT]. A recent 
experiment showed that dynamics of viral packaging 
is ultraslow and glassy resulting in significant hetero¬ 
geneity in packaging rates that vary from one virus 
to another m- It is noteworthy that the size of vi¬ 
ral DNA {N = (30 X 10^)/0.34 = 8.8 x 10^ bp) is 
only ^10 % of bacterial genome. Thus, the equilibra¬ 
tion time of DNA conformation based on reptation (or 
scaling) should occur 10 ^ times faster than in bacterial 
genome, which would contradict experiments [19]. To 
explain the ultraslow and heterogeneous dynamics of 
viral DNA packaging it is essential to consider the ef¬ 
fects of confinement, and our theory provides a natural 
explanation of the observations. 

Our study provides a general framework to quantify 
glassy dynamics of a polymer chain (a simple model 
for chromosome organization) and highlights the non¬ 
equilibrium aspect of a single polymer under strong 
confinement with clear implications for the variations 
in genome folding across different species. Dynami¬ 
cal implication of our finding 0bac ^yeast < ^ 

Shuman < ^virus, and thc Correlations of (/)yeast < 
with P{s) ^ for budding yeast [37] and < 

Shuman with P{s) ^ s~^ for mature human cells [9] 
provide a new framework for understanding the ori¬ 
gin of qualitatively distinct chromosome organization 
in various organisms and cell types. 

Given that cellular environment is replete with 
crowding particles, the volume fractions estimated 
here for different organisms may well be only lower 
bounds, and thus we expect that glassy dynamics is 
prevalent especially in higher-order organisms. To 
overcome topological constraints, fluidization or equi¬ 
libration of nuclear environment using topoisomerase 
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or metabolic activity would be sometimes necessary 
for biological systems to execute their functions m- 
Furthermore, it is noteworthy that although there is 
not significant difference in genome volume fraction 
between human embryonic stem cell (hESC) and 
mature cell [43], these two cells have distinct P{s) 
{P{s) for hESC, P{s) ^ for mature cell) 

(sun], which may be linked to substantial variations 
in metabolic activity or specific interactions with 
nuclear envelope depending on the cell maturity. 
Although our conclusions here do not consider the 
role of active mechanisms on genome organization, it 
is plausible that equilibration machineries exploiting 
active forces are required when chromosome dynamics 
is intrinsically glassy, as appears to be the case in 
higher organisms. 
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Supplemental Material 


Model. In order to assess the conditions describ¬ 
ing the onset of glassy dynamics of a confined flexible 
polymer we introduce a model in which the potential 
energy is given by 


N-l N-1 N N 

U{n, ...,nv) = +5] + E 

i=l i=l j=i-\-l i=l 



where a is a bond length, k = 80 /cpT, e = 2 /cpT, 
fi is the position of the monomer. To model the 
effect of confinement we placed the polymer chain in 
a sphere surface of radius Rg. The interaction be¬ 
tween the monomers and the sphere surface is repul¬ 
sive, given by the last term in Eq.(l). 

We performed Brownian dynamics simulations of 
a self-avoiding polymer under spherical confinement 
with N = 50, 75, 100, 125, 150, 200, 250, and 300 by 
integrating the following equations of motion, 

uV ' —> 

= -^r.U{n, ■ ■ ■ ,nv) + ri(i), (5) 


where Ti{t) is the Gaussian random force satisfying 
the fluctuation-dissipation theorem, (r^(t) • = 

6kBT(S{t — t')5ij. With the Brownian time defined as 
r = o?/D where D = ksTj^^ we chose the integration 
time step 6t = 8.6 x 10“^ r as a compromise between 
accuracy and computational cost. 

For eukaryotic genome, D = QTuha/ 2 ) ~ 

GttX (0.89X 10-^V/m2.sec) X lOnm ^ [IW? j S with T] = 

0.89 X 10“^Pa • s and a ^ 20 nm; and hence we set 
r = /D ^ 1.6 /as and 6t ^ 13.7 ps. 



FIG. SI: The protocol used to reduce the confinement size 
(Rs) for different N. The reduction rate of confinement 
size near the dynamical arrest point is similar for all N as 
r ~ —0.15(Rs/tt)/step. 


We gradually reduced Rg from ARg to the value at 
which Ta or Eqs. (2) and (3) in the main 

text) starts to diverge. Although the detailed proce¬ 
dure of reducing the confinement size Rg varies with 
N, the rate of Rg reduction r = ARg/At ^ 0.15 
(Rg/a) /(2 X 10^ X 5t) is almost identical for all N when 
Rg approaches to the point of dynamical arrest. The 
Rg values varied in the simulations are listed in Ta¬ 
ble |Ij and the time-dependent protocol of reducing Rg 
is plotted in Fig. SI. At each Ri^\ we simulated for 
2 X 10^ St and took the last conformation from the pre¬ 
vious simulation at R^^ as the initial conformation 
for simulation in Ri^\ We reduced Rg from R^^ to 
linearly for 2 x lO^St^ allocated the next 2 x 10^ St 
for an equilibration, and used the rest of 2 x 10^ steps 
to calculate Fq^^^{t) and We generated 10 inde¬ 

pendent trajectories for N < 150 and 25 for N > 200 
to improve the quality of statistics. 

It is worth emphasizing that the critical volume 
fraction (fc is robust and insensitive to the range of 
confining speed. To show this, we used two different 
confining speeds for a polymer with N = 150: one 
is r/ = — 0.30(Ps/a)/step, 2 times faster than r and 
the other is — 0.10(P^/a)/step, 1.5 times slower than 
r. Both from using Fq{t) and we obtained 

<f>c = 0.383 for both quenched and annealed cases, 
which is in full agreement with the regular case (see 
FigjSl). 

Volume fraction of a confined polymer. When 
a polymer is confined to a sphere with radius P^, the 
size of the polymer P^ can be related to the radius of 
gyration for polymer in free space (Pp via the follow- 
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N 

50 

75 

100 

125 

150 

200 

250 

300 

i — 1 

Rs — 10a 

10a 

10a 

10a 

10.8a 

13.3a 

13.3a 

13.3a 

2 

9.2a 

9.2a 

9.2a 

9.2a 

9.2a 

11.7a 

11.7a 

11.7a 

3 

8.3a 

8.3a 

8.3a 

8.3a 

8.3a 

10.8a 

10.8a 

10.8a 

4 

7.5a 

7.5a 

7.5a 

7.5a 

7.5a 

9.2a 

9.2a 

9.2a 

5 

6.7a 

6.7a 

6.7a 

6.7a 

6.7a 

8.3a 

8.3a 

8.3a 

6 

5.8a 

5.8a 

6.2a 

6.2a 

6.2a 

7.5a 

7.5a 

7.5a 

7 

4.3a 

4.3a 

5.8a 

5.8a 

5.8a 

6.7a 

6.7a 

6.7a 

8 

4a 

4a 

5.5a 

5.5a 

5.5a 

6.3a 

6.3a 

6.3a 

9 

3.7a 

3.7a 

5a 

5a 

5a 

6a 

6a 

6a 

10 

3.3a 

3.3a 

4.7a 

4.7a 

4.7a 

5.7a 

5.7a 

5.7a 

11 

3a 

3a 

4.3a 

4.3a 

4.3a 

5.3a 

5.3a 

5.3a 

12 

2.7a 

2.7a 

4a 

4a 

4.2a 

5a 

5a 

5a 

13 

2.5a 

2.5a 

3.7a 

3.7a 

3.8a 

4.7a 

4.7a 

4.7a 

14 

2.3a 

2.3a 

3.3a 

3.3a 

3.3a 

4.3a 

4.3a 

4.3a 

15 

2.2a 

- 

3a 

3a 

3.2a 

4.2a 

4.2a 

4.2a 

16 

2a 

- 

2.7a 

2.8a 

3a 

4a 

4a 

4a 

17 

- 

- 

- 

2.5a 

- 

3.7a 

3.7a 

3.8a 

18 

- 

- 

- 

- 

- 

- 

- 

3.7a 


TABLE I: Rs values used for the simulations with vari¬ 
ous N. In the table, a is the monomer-monomer distance. 
Rs was sequentially reduced from z = 1 to z = Zmax ac¬ 
cording to the procedure described in the SI text. The 
initial conformation was taken from the last conformation 
of the previous run except z = 1 where we generated un¬ 
constrained chain conformation. We set R^J^^ to four times 
the Rg of the unconstrained chain, and slowly decreased 
Rf'> to 




FIG. S2: Robustness of 4)c value for different reduction 
rate of conhnement size. In the legend, “Fast” denotes the 
faster conhning speed; “Slow” is for the slower one; and 
“Regular” is the speed used in Table 1. Black lines are 
the fits for the values obtained at the regular speed. In all 
cases, we obtained 0c = 0.383 for N — 150. 



FIG. S3: Mean spatial distance of polymer with N — 300 
as a function of intersegmental separation s for varying 
volume fraction 0. Log-log plot is shown on the right 
panel with the dotted line expected for the scaling of SAW 
{R{s) ~ Note that the condition of conhnement 

{Rg > Rs) trivially gives rise to the plateauing of R{s) 

[SlIIS]. 



FIG. S4: 0 value from the ht of ~ for 

the conhned polymer {N = 300) as a function of 0. It is 
noteworthy that the decrease of 0, the phenomenological 
stretching exponent that characterizes the extent of glassi¬ 
ness, is consistent with our observation that the polymer 
dynamics becomes more glassy with increasing 0. 


f{x) ^ x^. From ^ / Rg)^ ^ the exponent 

p ought to be p = {diy)~^ — 1. Therefore, substitut¬ 
ing Rg = aN^ where a is the Kuhn length, one gets 
R^g = Rs{a/Rs)^^^^N^^^. 

A definition of polymer volume fraction (0) using 
the ratio between Rg and 0 = (Rg/Rg)^ gives 
distinct scaling of 0 with N, depending on the strength 
of confinement: 


ing scaling relation with x = Rg/Rgi 

Rl = Rp{x). (6) 

(i) Under weak confinement {x ^ 1), correspond¬ 
ing to large Rg, the chain statistics will be unaltered 
Rgr-^Rgr^ with u = 3/5, and thus f{x) ^ con¬ 
stant. (ii) In contrast, a strong confinement {x 1) 
induces polymer collapse, so that Rg ^ and 


^ - (ky -1 

I {i^y^'' (strong, i?° » R,) 

(7) 

where Ijv = d for the case of strong confinement. 
Note that this definition of 0 is invariant under 
coarse-graining. 





















Radial distribution function. We used 

AT-l N 

= N(N-l) ^ ^ 

^ ^ i = l j = i + l 

to capture the extent of packing between monomers 
in Fig.lc. 

Contact probability. Contact probability as a 
function of genomic separation \i — j\ = 5 in Fig.Id 
is given by, 

Ei=i -5) 

(9) 

where ©(...) is the Heaviside step function. 0(x) = 1 
for X > 0; otherwise Q{x) = 0. 


statistics should obey that of self-avoiding walk. 
Given the distance distribution Ps{r) between two 
interior points separated by s along the contour, the 
contact probability is defined as P{s){^ Ps{r = 0)). 
From Ps{r) ^ {l/s^Yf{r/s^) ^ {l/s^Y{r/s^y for 
r <C s, where g is the correlation hole exponent 
and g = O 2 for two interior points [29]. The scaling 
exponent should be similar to the probability of two 
interior points of a SAW chain to be in contact, 
P{s) - ^ with d = 3, 6>2 = 0.71, 

u = 0.588 [29]43T] . In accord with this expectation, 
our simulation shows a = 2.18 in the absence of 
confinement {Rs/a 00). Note that for Gaussian 
chain (or polymer melt) ^ = 0, = 1/2, and d = 3, so 

that we retrieve the scaling relation for an equilibrium 
globule P{s) ^ in the above. 



e/ksT 


FIG. S5: (a) Snapshots of polymer under strong confine¬ 
ment (0 = 0.402). Monomers, colored based on the energy 
value, underscore the spatial heterogeneity of stress in the 
organization of the polymer, (b) Monomer energy distribu¬ 
tion, P{s), at 0 = 0.402 for different range of r: r/Rs < 0.5 
for core and 0.5 < r/Rs < 1 for the surface. Together with 
the snapshot displaying the interior of the globule on the 
right, P{£) for the different range of r highlights that the 
spatial heterogeneity of the monomer energy is present in 
the interior as well as on the surface of globule. 

Scaling relationship of contact probability 

for SAW. In the absence of confinement, the chain 



0 = 0.313 


# 

).402 ^ 

@ • 

FIG. S6: Formation of fractal-like globules from self¬ 
avoiding chain with increasing extent of confinement (0 = 
0 ^ 0.402). At 0 = 0.402, the globules display segregated 
domains with ultra-slow mobility. 






